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FOREWORD 

Astro Aerospace Corporation is conducting studies on design concepts for 
large space structures. This is one of a series of reports Issued as part of 
those studies under Contract No. NAS1-17536. 
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SECTION 1 
INTRODUCTION 


Meeting the power requirements of future space missions such as the Space 
Station will require advances in space power technology. Most currently oper- 
ating satellites have relatively low power requirements and use a low effi- 
ciency photovoltaic solar-array system to supply all of the power needs. 
Batteries are used during the part of the orbit in the earth's shadow. 
However, preliminary estimates of the power requirements for the proposed 
Space Station indicate that a photovoltaic power system would require an area 
of solar cells of thousands of square meters. Such arrays would have a large 
aerodynamic drag and would consequently affect the station operation and the 
selection of the station architecture. Station operation could be greatly 
enhanced and constraints on station architecture minimized by reducing the 
need for such large areas. Preliminary results from recent studies have shown 
the potential for significant area reductions using a solar thermal approach. 
In this approach, concentrated solar energy is used to heat a fluid which then 
drives a power-producing heat engine. 

The potential reduction in area using a solar thermal approach comes from 
two sources. One is the higher conversion efficiency of a solar thermal 
system. This advantage has been known for years and is controlled by funda- 
mental physical laws. The second Is In the overall power system operation and 
design. Recent studies have shown that the overall system operation would 
benefit by storing energy (for use while in the earth's shadow) as thermal 
energy rather than as electrical energy. The Increase In efficiency due to 
the combination of storage and conversion can reduce the surface area required 
for the Space Station power system by a factor of five. 

The advantage of the photovoltaic system is Its simplicity and 
insensitivity to geometric errors In orientation of the solar array. In 
addition, a large amount of development work has resulted from the frequent 
use of the technique. The structural needs and approaches have been studied 
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In detail and knowledge is well advanced on the weight, cost and performance 
that can be expected. 

In the solar-thermal approach, the solar radiation must be concentrated at 
a receiver so that high temperatures are produced. The amount of concentra- 
tion required is dependent on the type of power conversion chosen, but ranges 
from a value of 1,000 to 2,000. Such high concentration ratios place much 
more stringent demands on the concentrator surface and its supporting struc- 
ture than exist for the photovoltaic structural components. Consequently, the 
performance of a solar-thermal power system is very dependent on the perfor- 
mance and weight of the concentrator structure. 

During the 1960s, a large number of investigations were conducted on the 
solar concentrator and its support structure. This early work is reported in 
a variety of forms. Examples can be found in the Power Systems Conference 
Proceedings (References 1 and 2). Various types of reflectors are discussed 
in the next section. The type that received the most advanced development 
consists of honeycomb-core-stiffened reflector panels. A deployable collector 
with such panels as petals reached the stage of prototype hardware 
(Reference 3) . 

The success of photovoltaic systems for low power requirements, the lack 
of space missions requiring large power levels, and budgetary constraints led 
to the cessation of work on solar-dynamic (as well as nuclear-dynamic) power 
systems in the late 1960s. Not until the middle 1970s did interest rekindle, 
this time for terrestrial use, because of the disruptions in oil supply and 
price. In order to provide an alternative energy source, a multi-year program 
was undertaken, aimed at the construction of long-life, cost-effective 
paraboloidal-collector power-generation units. Several collector concepts 
have been developed and tested; the most successful consist of inexpensive, 
stiff reflector panels mounted on a truss supporting structure. (See 
Reference 4.) The current Space Station studies are benefiting from the 
technology developed for terrestrial application. (See Reference 5.) 

In the present paper, various types of solar concentrators are discussed 
from a structural point of view. The discussion includes reflecting surfaces 
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of a membrane type as well as stiffened panels. The methods for supporting 
the reflecting surface are described. The deep truss supporting structure Is 
identified as providing high accuracy at reasonable cost, although on-orbit 
assembly is required to mount the reflector panels to the truss. 

Ray-tracing analyses are performed to evaluate the tolerable amount of 
slope inaccuracy of the nominally paraboloidal reflecting surface of a 
concentrator. For the high concentration ratios (approximately 2,000) 
required for a Brayton-cycle heat engine, an rms slope error of three 
milliradians is acceptable. Next, a special nonparaboloidal reflector 
concept, composed of Identical square, spherically curved panels efficiently 
located, is analyzed by ray tracing. Such a configuration should reduce the 
costs of panel fabrication. 

Because the identical-panel concept with large panel sizes is shown to be 
feasible, a configuration is generated which combines the square panels with a 
Pactruss supporting structure. The resulting design Is described and Its mass 
breakdown is estimated. 
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SECTION 2 

SOLAR CONCENTRATOR CONCEPTS 

In this section, attention is devoted to those concentrator configurations 
which are so large that they must be carried Into orbit in a packaged 
condition. The collecting area is greater than about 100 square meters. 
Therefore, fixed dishes, which can be assembled and aligned on the ground 
prior to launch, are precluded. 

The concentrator must not only be packaged, but also be compactly stowed 
for launch. The Space Station will require numerous Shuttle launches; the 
concentrator must not add unnecessarily to the number. 

In general, a solar concentrator can be considered to be composed of three 
components: the collecting surface, the supporting structure, and the mounts 
which connect the concentrator to the rest of the spacecraft. It is conven- 
ient to consider each component separately. 

Collecting Surface 

Lenses - The surface can consist of a lens or a reflector. Lenses can be 
constructed of thin layers by using the Fresnel principle. Indeed, the lens 
can be structurally a flat membrane and supported as a drumhead. The lenses 
have the advantage of being relatively insensitive to geometric errors. On 
the other hand, the energy must be transmitted through the lens material. 
Maintaining low loss throughout the long-term exposure to the space environ- 
ment places stringent requirements on the lens material. In addition, the 
supporting structure will block the energy unless great care is taken. 

Reflectors - The reflecting surface has been adopted almost universally 
for solar concentrator for terrestrial and space applications. The surface 
must be a highly specular reflector and must be oriented accurately in order 
to direct the reflected rays into the thermal receiver. Analyses such as 
those described in Section 3 show that the orientation accuracy required is 
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not as great as that for many of the high-performance RF missions presently 
being undertaken. The sun subtends about nine milliradians in arc. Little is 
gained by demanding reflected-ray accuracies better than about three 
mi 111 radians, even for a concentration ratio of 2,000. 

The achievement of highly specular, low-loss reflection for long-term 
space missions is a more difficult objective. Even with solid surfaces, about 
ten percent of the incident energy is lost; with metallized films the specular 
reflectivity is reduced to perhaps 80 percent after long exposure. The 
reflection losses are minimized in most collector arrangements by avoiding 
multiple reflections. 

Film Surfaces - Film-type reflector surfaces are attractive from packing 
and weight standpoints. Furthermore, only the film surface affords the possi- 
bility of uncomplicated deployment without requiring on-orbit assembly. 
However, the surface must be doubly curved with positive Gaussian curvature. 
In order to achieve this, some method must be found to apply a pressure across 
the membrane. Inflation will accomplish this objective, but leakage will 
require prohibitive amounts of make-up inflatant unless the collector is very 
large (kilometer scale), the pressure is low, or the pressurization time is 
small. The last alternative is probably the most feasible. Inflation is used 
to stretch the wrinkles from the packaged reflector. The resulting reflector 
is then stiff enough to maintain Its shape, either because the surface mate- 
rial is a multilayered film or because it is a composite which is cured on 
orbit after inflation. 

In any event, the film will operate at very low stress levels. Thus, the 
shape accuracy must be attained primarily by enforcing the correct geometry 
during fabrication. As is shown In the next section, efficient solar concen- 
trators tend to be deep dishes; for these the in-surface stiffness dominates 
the structural behavior. 

An alternate method of applying the shaping pressure is by electrostatic 
attraction. This technique avoids the gas leakage problem. It furthermore 
offers the opportunity of controlling the loading and hence, the shape. The 
latter controllability can be useful for shallow dishes, but not for the types 
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needed for solar concentrators, unless the surface Is divided Into numerous 
panels. 

Panel Surfaces - The alternative to the membrane reflector surface Is an 
assemblage of stiff panels that are Individually small enough to fit In the 
launch vehicle. Each panel Is considered to be stiff enough to maintain Its 
own shape. 

There are two primary textures for the panel construction. One is the 
honeycomb sandwich; the other Is the thin shell, stiffened or monocoque. If 
the shell Is monocoque, the panel boundary must be dimensionally stable in 
order to preclude warping. 

Panels can be constructed from a wide range of materials. Long-life sta- 
bility of the reflective and dimensional properties would appear to favor 
metallic construction. The Sunflower petals (Reference 3) were essentially of 
aluminum construction, face sheets and honeycomb core. Recent work has dealt 
with graphite- and Kevlar-reinforced composites for the face sheets. 
Producing a smooth reflecting surface will probably require a shaping mandrel 
in either case. 

Masses - The film surface has a unit mass of the order of 0.1 kg/m 2 , 
whereas the panel concept will weigh at least one kg/m 2 . If very precise 
surfaces are required, the mass can Increase by another order of magnitude. 

Supporting Structure 

Film Surfaces - Concentrator configurations utilizing film reflector sur- 
faces usually consist of a transparent pressure containment surface joined to 
the reflector at the rim. The containment surface can be shaped symmetrically 
to the reflector so that the pressurized vessel Is lenticular. [See Figure 
1(a) and (b) taken from Reference 6.] The rim must be capable of carrying 
compression without local crippling. It must also be stiff and stable enough 
to maintain the shape of the edge after the internal pressure Is removed. 
Rims for inflated reflectors tend also to be inflated. Deployment is achieved 
by pressurizing in a controlled fashion. 
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If electrostatic forces are used to tension the membrane, no frontal mem- 
brane Is needed. On the other hand, a substantial back-up truss Is required 
to support the electrostatic devices. In addition, the membrane and the 
charged devices must be shielded from the space plasma In order to avoid elec- 
trical leakage. 

Panels - The supporting structure for panels can be made up In a variety 
of ways. At one extreme, the panels are stiff enough to serve as their own 
structure. At the other, a separate back-up structure (a truss, for example) 
is furnished and the panels are attached to it. An Intermediate configuration 
could be one In which the panels derive their main support from their own 
stiffness but are externally supported at several points, say at the rim, in 
order to remove large-scale errors. 

The aforementioned Sunflower is an example of the self-supporting panel 
approach. Each petal Is hinged at its inner edge to a central hub. When 
deployed, the petal Is cantilevered at Its inner edge and Is Independent of 
its neighbors. This version Is about 30 feet in diameter. For larger 
collectors, more folds become necessary and adjacent panels must be 
structurally joined. Some possible concepts are discussed In Reference 7. An 
example is shown as Figure 2. Of course, the objective of these designs Is to 
enable automated deployment. The more In-space attachments that are required, 
the greater is the complexity of the deployment apparatus. 

Deployable configurations using nearly hexagonal panels are possible as 
discussed, for example, in Reference 8. By using some secondary folds, a 
seven-hexagon collector about 12 meters in diameter can be stowed in the shut- 
tle cargo bay. If an Aft Cargo Carrier were available, almost 20 meters in 
aperture would be possible. (See Reference 9.) 

While many practitioners are convinced that the self-supported-panel 
approach can yield the accuracies necessary for even high concentration ratios 
and large sizes, it is the conviction of the authors, based on flight hardware 
experience, that the unit cost of this technique gets prohibitively large for 
large sizes. The ideas are illustrated qualitatively in Figure 3. The unit 
cost of the self-supporting-panel approach increases as collector size 
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Increases because of the added complexity of the deployer and, more 
importantly, the cost of the high precision necessary to maintain required 
accuracy. In addition, the cost of testing and proving the accuracy of the 
system on the ground becomes high. There is, in fact, a practical limit on 
the possible size of this approach. 

On the other hand, the truss-supported approach costs more for small 
aperture sizes because of its Increased complexity and the costs of assembling 
in orbit. But when the size Increases, the unit cost decreases because much 
of the additional cost of the truss support is not very size dependent. The 
inherently stlffer structure with greater depth Is less sensitive to fabrica- 
tion errors and can be tested and measured on the ground more easily. The 
result is that much larger apertures are made practical by the truss-support 
approach. In addition, that approach could be less expensive even In the 
smaller size range considered in this paper, if high concentration ratios 
(approximately 2,000) are required. This may be particularly valid for the 
Space Station If manned or robotic assembly becomes a routine matter so that 
the costs of such assembly were reduced. 

Examples of the truss-support concept are shown In Figures 4 through 7 
taken from Reference 10, Figure 8 from Reference 11 and Figures 9 through 11 
from Reference 9. They Include a large range of implementation approaches 
from automated deployment, through assembly with a remote manipulator, to 
extravehicular assembly of Individual panels to the support truss. Most are 
unproven but all are reasonable and obey the laws of physics. 

The approach in Figure 4 Is to attach panels to a previously deployed sup- 
port truss. In Figure 5, the truss Is attached to the panel before launch. 
The units are deployed and assembled to their neighbors on orbit. 

An automated deployment scheme Is shown in Figure 6 and 7. In this 
approach, an intelligent robotic deployer expands each truss module and 
attaches it to its neighbors as needed. The modules are hinged together so 
that the number of attachments is minimized and the hinges help to guide the 
automated construction. 

A very simple deployable concept is shown in Figure 8. This approach is a 
descendant of the Extendible Support Structure which was flight-proven on the 
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Synthetic Aperture Radar on the SEASAT spacecraft. Note that the panels nest 
when packaged so that deep dishes can be stowed efficiently. 

Figures 9, 10 and 11 show an assembly-intensive approach where large panel 
modules with pre-attached stowed support-truss segments are carried to orbit 
in the After Cargo Carrier. Each module Is removed from the ACC, deployed and 
assembled. Collector areas of 300 square meters are possible with this 
approach. 

Masses - The mass of the support structure can vary from zero with the 
self-stiffened panels to as much as the panel mass. For film reflectors, the 
support structure mass is equal to that of the reflector if a frontal membrane 
is used and will be many times greater (perhaps ten kg/m 2 ) if electrostatic 
pressure is used. 


Mounting Structure 

The collector must be positioned and oriented accurately with respect to 
the receiver which gathers the concentrated solar power. Furthermore, the 
combination must be pointed at the sun to an accuracy of a couple of 
milllradlans. And the needed structure must not block too much solar energy. 

For the inflatable collector, It is tempting to connect the frontal mem- 
brane to the receiver. This can lead to great difficulty with handling the 
heat which is re-emitted from the receiver cavity. Furthermore, such a con- 
figuration would require continuous pressurization. Hence, the usual concepts 
provide struts, perhaps inflated, connecting the receiver to the rim. Hard 
points for connections and means to distribute the point loads into the col- 
lector must be provided. 

The self-supported panel concepts also are subject to the difficulty of 
providing hard points for attachment. On the other hand, attachment to the 
truss support structure is straightforward. Note, however, that space-use 
trusses are generally built of slender members so that care must be taken to 
ensure that the attachments cause no local moments to be applied to the 
joints. 
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The attachment -struts can be long, hinged members (see Figure 7, for 
example) or deployable or erectable beams. The choice Is usually dictated by 
the texture of the collector, the type of spacecraft Interface, and the stow- 
age space available. The designer should avoid placing the struts close to 
the mouth of the receiver lest the heat damage the structure. 

Selected Approach 

The approach selected for further study herein Is the truss-supported 
panel concept. The main reason Is Its potential for constructing large col- 
lectors with high concentration ratios. Secondly, the Space Station is 
assumed to be able to provide on-orbit assembly in a routine fashion. 

The cost-saving potential for the truss-supported approach has been 
pointed out. In addition, the structural integrity does not depend on the 
panels themselves. This not only allows more latitude in panel design but 
also simplifies the task of replacing panels on orbit. The higher natural 
vibration frequencies enabled by the deep truss also simplify operational 
procedures. 

In the past several years, a great amount of effort has been devoted to 
the large infrared telescope which has an aperture of 20 meters. One of the 
attractive configurations is described in Reference 12 and is shown in Figure 
12. The panels (in this case, weighing around ten kg/m 2 ) are mounted on a 
newly-invented synchronously-deployable truss structure called Pactruss. The 
deployment kinematics shown in Figure 13 have the property that the bays of 
the truss are very strongly coupled together. Thus, reliable synchronous 
deployment can be achieved with only a few actuators. A photograph of an 
early working model is shown in Figure 14. 

The Pactruss deployable structure can also be configured in a square 
pattern. Such an example, in which the reflecting surface is offset, is shown 
in Figure 15. Reference 11 contains more detail on the square-cell geometry. 
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SECTION 3 

EFFECTS OF GEOMETRIC ERRORS ON CONCENTRATOR EFFICIENCY 

Collectors of large effective concentration ratio are sensitive to small 
errors In the shape of the reflecting surface and to small errors In pointing 
the axis of the concentrator at the sun. The finite size of the sun's image 
also contributes to a reduction In the effective concentration ratio. Errors 
of a different sort arise when the reflecting surface consists of a mosaic of 
small independent spherical panels mounted on a common supporting truss 
structure. The effects of such errors on the effective concentration ratio of 
the concentrator are reported in this section. Note that the surface Is 
assumed locally to be a perfect specular reflector in this section. 

Effects of Random Slope Errors 
on the Effective Concentration Ratio 
of a Contlnous-Parabololdal Reflector 

The effect of geometric imperfections In the reflector surface on concen- 
trator efficiency is a topic which has been studied by many Investigators 
(see, for example, Reference 13). The objective of the study presented herein 
is to focus on the effects of those Imperfections which are of particular 
Importance In the structural design of precision reflector surfaces and sup- 
porting structures. 

Consider first a perfect paraboloidal reflector with axlz z, as shown in 
Figure 16. The reflecting surface Is described by the equation 

r 2 

z * — (1) 

4F 

where r is the distance from the z-axls to a point on the surface, and F is 
the focal length, or distance from the origin 0 to the center of the 
receiver R. For purposes of analysis, the receiver Is assumed to be a circu- 
lar disk of diameter d and oriented orthogonal to the z-axls as shown In the 


11 


figure. The diameter of the frontal projection of the paraboloidal reflector 
is D. 


For a reflecting surface defined In this way, any Incident ray parallel to 
the z-axis and Intersecting the reflecting surface at a point q will reflect 
into a ray which Intersects the receiver at point R. However, the rays from 
the sun incident at q are not all parallel to the z-axis since the sun's image 
is a disk in the sky subtending an angle 2e $ . At a distance of approximately 
one astronomical unit, the sun half-angle e $ Is 

e $ = 4.5 milliradlans (2) 

Thus, as reflected at a point q, the rays from the sun fill a right circular 
cone with apex angle 2e and axis qR. 

In order to capture all of the sun's rays reflected at q, the receiver 
diameter d must be selected to be sufficiently large. However, d should not 
be chosen any larger than necessary to capture all the sun's rays emanating 
from every point q on the reflector of diameter D. To do so would reduce the 
geometric concentration ratio, defined as C where 

/ D\ 2 Frontal Area of Reflector 

C h - — (3) 

\ 2 / Area of Receiver 

and would allow needless power losses in radiation from the hot receiver. The 
task of selecting the proper d for a given D (or vice versa) requires consid- 
eration of the geometry of the Intersection of the cone of reflected rays and 
the receiver disk. For this purpose, consider a planar coordinate system cen- 
tered at R and containing the receiver disk, as shown in Figure 17. In this 
coordinate system, the receiver is a circle of diameter d centered at the 
origin. 

The image of sun's rays on the receiver plane is an ellipse. Let p repre- 
sent the distance from R to q. Then the major and minor axes of the ellipse 
are 2pe /cos 0 and 2pe , respectively, where 0 is the angle between Rq and RO, 

3 J 

as shown in Figure 16. Let the u and v axes be defined as parallel to the 
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minor and major axes of the ellipse, respectively, as shown In Figure 17. For 
a reflector surface without geometrical Imperfections, the ellipse is centered 
at R. However, in the presence of small errors in the slope of the reflector 
surface at q, the ellipse Is not concentric with the receiver disk. In that 
case, the center of the ellipse is translated a distance u„ and v 0 where 

u„ = 2pa u cos 0 (4) 

v 0 = 2pa y /cos 0 (5) 

where a u and a v are the slope errors at q in the u and v directions, 

respectively. 

If the sun's image lies entirely within the receiver disk, then all of the 
incident energy at q is delivered to the receiver. However, in the more gen- 
eral case, a portion of the sun's image lies outside the receiver disk as 
shown in Figure 17, and only a fraction of the incident energy at q is deliv- 
ered to the receiver. If the sun's image on the plane of the receiver is 

assumed to have uniform density, this fraction is the ratio of the area of 
intersection of the ellipse and circle in Figure 17, to the area of the 
ellipse, and may be called the local "capture ratio" y at q. The spatial 

average y of the local capture ratio over the projected frontal area of the 

reflector then may be used to define the "effective concentration ratio" 
in terms of C as follows 

C,ff - VC (6) 

Note that C„.r* can be considered to be the number of "suns" concentrated at 

eff 

the receiver. It is generally smaller than the geometric concentration ratio 
C defined in equation (3). It is which is of primary importance in solar 
energy collection. 

To investigate the effects on of random errors in the slope of the 
reflector surface, consider the special case when a u and a y are Independent 
Gaussian random variables with zero mean and identical standard deviation S Q . 
In this case y is also a random variable, and so is C gff , although C is 
deterministic. The expected value E [C gff ] of C gff is then determined by 


13 



(7) 


E [C eff ] = CE[Y] 

In order to determine the expected value. It Is necessary first to calculate y 
as a function of a u , a v , 0, F/D, and C. Numerical evaluation of y by simple 
strip-theory quadrature was used to develop the results presented herein. 
Averaging y over the reflector, again by numerical quadrature, yields Y. 

For the situation just described, numerical results for E[y] were obtained 
by an additional numerical quadrature of the equation 


E[Y] 


F 1 

Y (a , a , C) exp. 

U V D 2tt $ 2 



+ O 


2 S 0 2 


da u da v 


( 8 ) 


To facilitate the quadrature, the domain of integration in equation (8) was 
mapped to the strip p > 0, 0 < 4» < u/2 by utilizing symmetry and with the sub- 
stitution 


a u = /2p~ cos <j> 

(9) 

a y = /2p~ sin <J> 

(10) 


Note that the domain becomes finite since y is zero for large enough values of 
p. Results are shown in Figures 18 through 22 as plots of C versus E[C eff ] 
for different rms slope errors o Q , each figure corresponding to a different 
value of F/D. 

The effects of large o Q are typified by the curves in Figure 18. For 
example, when o Q = 10 milliradians, the results show that for F/D = 0.3, it is 
not possible to obtain a mean effective concentration ratio larger than about 
1700, no matter how large a geometric concentration ratio C is chosen. 
Furthermore, the same figure shows that an E[C gff ] of 2,000, for instance, 
requires a C of about 3,000 when o Q is 5 milliradians but only a C of about 
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2,250 when o Q Is 3 mi Hi radians. For the range of EC c e ff] shown in the 
figure, the curves for o Q ■ 0 and o Q = 1 milliradian are indistinguishable. 
For E[C ff ] up to 2,000, the values of geometric concentration ratio are 
essentially equal to E C c e ff3* Similar basic trends for other values of F/D 
ranging from 0.4 to 0.7 are shown in Figures 19 through 22. 

From the design viewpoint, it is most often the case that a predetermined 
value of is required. Then the design is specified in an overall geomet- 
ric sense by selecting the shape F/D and geometric concentration ratio C. It 
is therefore useful to consider the results in the form of C versus F/D for 
various levels of surface error o Q , as shown in Figures 23 through 25. Each 
figure corresponds to a different level of E[C gff ]. These figures again show 
that a smaller C is required for smaller o Q than for larger o Q . However, the 
figures also show the sensitivity of the required C to F/D. For example, when 
E[C e ff] is required to be 2,000, Figure 25 shows that an optimal choice of F/D 
exists. When o Q is 5 milliradians, the curve displays a minimum C at an F/D 
of about 0.45. For values of F/D substantially different from this optimal 
value, a significant increase in C results. However, for smaller values of 
rms surface error, the optimum value of F/D increases, although the curves 
show very little sensitivity to reasonably large deviations from the optimal 
F/D. 


For desired mean effective concentration ratios up to 2,000, the surface 
slope error can be three milliradians without noticeable penalty for F/D s 

0.5. 


Effects of Location, Orientation and Size of Nominally 
Square Spherical Panels on the Effective Concentration Ratio 

The application of large solar concentrators in space presents some diffi- 
culties in transportation. One practical approach to overcoming such diffi- 
culties is to erect the concentrator in space as a mosaic of panels attached 
to a supporting truss structure. For such a panel design, obvious simplifica- 
tions are gained by making the panels identical in shape, thus deviating from 
the desired paraboloidal surface. Near the center of a paraboloid, a very 
close approximation to the surface is provided by a sphere of radius 2F, where 
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F is the focal length of the paraboloid. However, at points far from the 
center, such spherical panels may not match the paraboloidal surface very 
well. In order to minimize the degradation in effective concentration ratio 
caused by the use of spherical panels, various strategies for the location and 
orientation of the panels may be considered. 

In order to effectively assembly a mosaic of panels without either exces- 
sive overlap or large holes, consider each panel to have a plan form which is 
a square of side length "a," as shown in Figure 26. Furthermore, let the 
radius of curvature of the panel be 2F. 

Consider an x, y, z coordinate system with the z-axis corresponding to the 
reflector axis, and with origin 0 at the center of the reflector surface, as 
shown in Figure 27. The center of the receiver R is located a distance F 
along the z-axis. 

Let the panel be oriented so that for small angles of tilt, its projection 
onto the x, y plane has sides which are always parallel to the x, y coordinate 
directions, as shown in the figure. Furthermore, let the panel center be 
located at c, and let the projection of c onto the x, y plane be c'. Then the 
points 0, R, c and c' lie in a plane containing meridional, lines of the 
reflecting surface, and labelled "meridional plane" in the figure. The 
meridional plane makes an angle with the x-axis. Let the angle between cR 
and the z-axis be 0 C , as shown in the figure. 

Consider a point p on the panel surface and also on the meridional plane. 
The location of p is determined such that an incident ray at p which is paral- 
lel to the z-axis reflects into a ray which intersects the receiver at its 
center, R. The distance from R to p is p p , and the angle between the ray pR 
and the z-axis is 0 p , as shown. In this case, the local normal to the panel 
surface at p makes an angle 8 p /2 with the z-axis. 

For given values of p p , 0 p and \|>, the center s of the generating sphere 
from which the panel surface is derived may be located from the equations 
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where x , y , and z are the coordinates of s. Note that s lies in the 
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meridional plane. The coordinates x q , y q , and z q of a generic point q on the 
spherical surface are then governed by the equation 



(12) 


For a given receiver diameter d, the location of the panel center c is 
determined in the following manner. First, consider temporarily locating the 
panel so that its center c coincides with p. While in this position, the 
panel size a is increased until the reflected rays from the worst case loca- 
tion on the panel (always a panel corner) contain sufficient error in direc- 
tion that they first intersect the perimeter of the receiver disk. Clearly, 
the reflected rays from any smaller panel would all lie entirely within the 
receiver disk, but for a larger panel, a fraction of the reflected rays would 
miss the receiver altogether. The value of this maximum panel size is noted. 

Next, the entire panel is rigidly rotated about the center s of the gener- 
ating sphere an incremental amount around an axis orthogonal to the meridional 
plane. During this motion, the panel sides remain parallel to the coordinate 
axes, but the center c is translated along a meridional line away from p, as 
shown in the figure. In this new position, a new maximum panel size is deter- 
mined by the procedure just described, and note is taken of the value. Then, 
the panel is similarly rotated through an additional incremental amount, and 
the process is repeated. The final location of the panel center is chosen as 
that location which permits the maximum possible panel size a m to be used. 
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It Is necessary at this point to discuss the nature of the errors Included 
in the analysis of the reflected rays. First, the effects of the finite size 
of the sun's image in the sky, as measured by e $ , are Included in a manner 
similar to that discussed in the previous section. In addition, a worst case 
mean pointing error in the alignment of the z-axis with the center of the sun 
is assumed. This angular pointing misalignment is e p in magnitude, and 
assumed for the Illustrations presented later to be 

e p = 0.15° = 2.62 milliradians (13) 

When analyzing the effects of e s and e p on the reflected rays at any point q 
on the surface, the local worst case for the direction of the pointing error 
is assumed. In addition to the errors just described, geometric errors are 
naturally introduced from the fact that the spherical reflecting surface devi- 
ates from the desired paraboloidal surface. However, the spherical surface is 
assumed to contain no random imperfections in this case. 

The three remaining parameters which affect the maximum allowable panel 
sizes are v|/, 0, and p p . For panels located on a nominally square grid aligned 
along the x, y coordinate directions, the panels will be located at various 
values of i|». The effects of on the maximum allowable panel size will be 
presented later in the results. 

For panels near the reflector axis, 6 p (and hence also 0 C ) will be small, 
but for panels near the outer rim, 0 p will be maximized for a given reflector. 
For spherical panels considered here, errors build up rapidly as 0 p is 
increased, and it is the panels near the rim of the reflector which govern the 
maximum panel size which may be used throughout. Consequently, only panels 
located near the rim are considered in the results presented later. 

The distance p p from the receiver center R to the point p on the panel 
surface has a major effect on the maximum allowable panel size. In general, 
one would expect the optimal choice of p p for a given 0 p would be that value 
which results in a best fit of the spherical panel surface to the local sur- 
face of the paraboloid at 0 p . However, selection of this best fit location is 
complicated by the fact that the spherical panel has a radius of curvature of 
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2F which is the same in all directions, whereas the paraboloidal surface has 
different radii of curvature in the meridional and circumferential directions. 

Three different approaches for selecting 0 p were investigated and are 
reported in the results. For a given 0 p , these are (1) selecting p p such that 
the radius of curvature of the spherical panel matches the local meridional 
radius of curvature on a paraboloid, (2) selecting p p such that the panel radius 

of curvature matches the local circumferential radius of curvature on a 
paraboloid, and (3) selecting p p as the geometric mean of the distances deter-* 
mined in (1) and (2) above. These three approaches result in the following 
selections for p n , labelled to correspond with the cases outlined above 

.... "Meridional Curvature" (14) 


. ... "Circumferential Curvature" (15) 


.... "Geometric Mean Curvature" (16) 

Results are presented for each of these three approaches for determining p p . 

Shown in Figure 28 are results for the dependence of the dimensionless 
maximum panel size (a ffl /F) as a function of the angular location 0 C of the cen- 
ter of a panel. Different curves are shown for different values of v|/. All 
curves in the figure are for the case where C g ^ = 2,000 and (p p /F) = 1. 
Comparing the curves in the figure, each displays a local peak of (a m /F) at an 
intermediate value of 0 . For very small 0 (corresponding to a high value of 
effective F/D), the value of (a m /F) in all cases is very small because of the 
•effects of finite sun size e . For large 0 (corresponding to small effective 
F/D), the value of (a m /F) in all cases is again small, this time because the 
spherical shape of the panel differs considerably from the shape of a 
paraboloid surface. 
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It Is further noted in Figure 28 that the case where ^ = 0 degrees (panel 
sides parallel to the coordinate axes x and y) is somewhat of a best case in 
that significantly larger panels may be used than for other values of ty. In 
many cases (although not in the case shown in the figure), the worst case 
occurs for ^ * 45 degrees. This worst case therefore governs the size of pan- 
els in a practical design, and remaining results are presented for v|> = 45 
degrees. 

The effects of pivoting the panel about the center s of the generating 
sphere so that c and p are coincident are shown in Figure 29. Shown in the 
figures are curves for (a m /F) versus 0 C for the case where (p p /F ) = 1 and 
C ef .f = 2,000. The solid line corresponds to results obtained by the pivoting 
maneuver previously described, while the dashed line indicates similar results 
obtained without pivoting so that p and c are always coincident. It is noted 
that the pivoting maneuver results in a significant increase in the maximum 
allowable panel size (a m /F). 

Presented in Figures 30 and 31 are results which show the effects of the 
three different strategies for locating the panel center. Note that in both 
figures (corresponding to = 1,000 and 2,000, respectively, with a worst 
case ^ = 45 degrees), the best choice is simply p p = F, which corresponds to 
locating the panel centers on a sphere centered at the receiver center R. The 
worst case in both figures is the choice (p p /F) = 1/cos (9 p /2) , which corres- 
ponds to locating the panel centers so that the panel radius of curvature 2F 
matches the circumferential radius of curvature of the local paraboloidal 
surface. Note further that in each figure there exists a local peak for the 
curve (a m /F) versus 0 C corresponding to the best choice p p = F. Hence, there 
exists an optimal effective F/D for each case. In the case where = 1,000 
shown in Figure 30, the peak occurs at 0 C a 1. In the case where C ef ^. = 2,000 
shown in Figure 31, the peak occurs at 0 C a 34 degrees, corresponding to 
F/D a 0.89. 

The foregoing results neglect the effects of surface slope errors in addi- 
tion to the pointing error of 0.15 degrees. Inasmuch as the analysis is based 
on the requirement that the entire reflected beam from all parts of the 
reflector is captured by the receiver, it is probable that the losses due to 
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one to three mllliradians of slope error would be small. The combined analy- 
sis is relegated to the future. 


Summary of Findings 

Based on the analyses presented in this section, the conclusion is that, 
for continuous paraboloidal reflectors with mean effective concentration ratio 
up to 2,000, an rms surface slope error as large as three milliradians may be 
tolerated without significant loss in concentrator efficiency. Furthermore, 
the optimal value of F/d in this case appears to be slightly less than 
0.5. 

For reflectors consisting of a mosaic of individual spherical panels of 
side length a m , the panel size may be as large as F/ 10 without significant 
loss in efficiency for mean effective concentration ratios up to 2,000. These 
conclusions are based on an analysis of the effects of errors due to aberra- 
tions in the mosaic reflector surface, finite sun size, and an overall con- 
centrator pointing error of 0.15 degrees. 

Furthermore, the results should not be very sensitive to moderate 
Increases in the pointing error due to local distortions, because the results 
are based on an assumption of complete capture in the worst-case direction for 
the worst-case corner on the worst-case panel in the mosaic. Therefore, only 
a moderate loss in captured energy would be incurred on a small portion of the 
surface area if a larger pointing error were present. 
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SECTION 4 

EXAMPLE OF TRUSS-SUPPORTED CONCEPTUAL DESIGN 


The analysis In Section 3 establishes that hlgh-efflclency, high- 
concentration ratio solar concentrators can be constructed of Identical panels 
that are properly placed and oriented. In particular, a concentration ratio 
of 2,000 is comfortably achieved with identical square panels having a size of 
one-tenth the concentrator diameter located at a distance from the receiver 
equal to the diameter. The conceptual design described in this section is 
based on those ratios. 

The choice of square facets for the Pactruss enables significant simplifi- 
cation in the design of the truss. Let the coordinate along the solar axis be 
z and the shape of the truss surface to which the panels are mounted be given 
by the form 

z = f(x) + f (y) 

The savings in joint similarity and tooling costs are obvious. 

Of course, the paraboloidal surface obeys the foregoing equation. A 
spherical shape does not, but the departures for the proportions considered 
herein are small. Furthermore, although the central points on the panels are 
supposed to be equidistant from the receiver, the outer corners, which are 
assumed to be mounted discretely at the truss nodes, do not lie on a sphere. 
A computer program (SOLARPAC.C) for calculating the function f(x) for chosen 
proportions of the panels and truss, the focal length, and the number of pan- 
els is included in the Appendix. This program considers panels oriented along 
the x-axis (or y-axis) whose central points are located on a sphere and whose 
tilt is adjusted so that the ray reflected from the central point is directed 
to the center of the receiver entrance. (The geometrical changes required to 
account for the improved accuracy discussed in Section 3 are small and are 
ignored here.) 
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In using SOLARPAC, the panel size must be chosen smaller than the strut 
length. Otherwise, Interference occurs between the corners of adjacent panels 
along the x = y direction. The resulting panel geometry Is shown In Figure 
32. The proportions are: 

Panel size = 2.0 meters 

Strut length = 2.15 meters 

Focal length = 21.5 meters 

Note the close meeting of the outer panel corners near x = y and the gaps 
elsewhere. 

Some consideration was given to the use of an offset concentrator geometry 
In order to avoid problems with blockage. But the loss of efficiency Involved 
in such an approach was considered to be unacceptable. This is the price to 
be paid for being able to use the Identical-panel concept. Note that the 
blockage problem Is avoided by providing a slot along the concentrator 
centerline in which mounting structure and radiator panels can be located. 
The truss structure carries through but the panels are omitted. 

The overall conceptual design as applied to the Space Station is shown in 
Figure 33. The concentrator is composed of 86 Identical two-meter-square 
panels, fabricated with a spherical radius of curvature of 43 meters. The 
capture area is 328 square meters. Each panel Is mounted directly to a 
Pactruss node at Its outermost corner; short standoffs connect it to the other 
nearby nodes. The support truss is attached to the receiver by means of an 
erectable mounting truss made of the same diameter and similar fittings as 
those of the primary Space Station truss. The mounting truss has the same 
five-meter dimensions In one plane but has a depth of only two meters In order 
to avoid blockage. The mounting truss Is connected to the Space Station Beta 
joint with a transition truss which has a variable geometry so that fine 
pointing is enabled. The attachment location is selected so that the center 
of gravity of the entire solar power generator Is In line with the transverse 
beam center. 

The radiator Is closely coupled to the receiver-converter in order to min- 
imize the length of heat lines. The radiator is also located in the slot and 
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is supported against out-of-plane motion by thermal-expansion-tolerant joints 
with the support truss. 

A possible scenario for constructing the new concept on the Space Station 
would be as follows: 

(a) Erect the mounting truss from the Beta joint. 

(b) Attach the stowed Pactruss to the mounting truss. 

(c) Deploy the Pactruss using its self-contained joint actuators. 

(d) Mount the reflector panels onto the Pactruss joints. 

(e) Attach the receiver/converter to the mounting truss. 

(f) Attach the radiator to the receiver/converter and the Pactruss. 

The resulting configuration is stiff and efficient. Its free-free funda- 
mental vibration frequency is estimated by the method in Reference 10 to be 
greater than ten Hertz. Its construction makes straightforward use of the 
assembly facilities which must be available on the Space Station. Estimates 
of the mass are as follows: 

Panel s : one-inch thick aluminum honeycomb with 

metallized graphite faces. 

344 m 2 @ 2.0 kg/m 2 688 kg 

Pactruss : one-inch diameter graphite tubes with 

0.020-inch wall. 2.0 joint factor. 

2,430 m @ 0.15 kg/m 355 kg 

Connecting 

Truss : two-inch diameter graphite tubes with 

0.060-inch wall. 1.5 joint factor. 

449 m @ 0.66 kg/m 296 kg 

The total mass of the collector and the mounting structure is 1,339 kilograms 
for a weight density of 4.1 kg/m 2 based on a frontal area of 328 square 
meters. 
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The stowage volume of the structural parts of the configuration is small. 
The volumes are estimated to be: 


Panels : 86 panels with three-centimeter depth 

per panel. 

2 m x 2 m x 2.6 m 

Pactruss : stowed size of 3-1/2 diameters per bay. 

1 m x 1 m x 4.3 m 


Connecting 

Truss : 100 struts of two-inch diameter, up to 

five long. 32 balls. 


0.5 m x 0.5 m x 5 m 


10.4 m 3 
4.3 m 3 


1.3 m 3 
17.0 m 3 
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SECTION 5 

CONCLUDING REMARKS 


The study described herein needs to be extended to Include more quantita- 
tive evaluation of the dependency between cost and required precision for the 
concentrator structure. The emphasis during this evaluation should be on a 
comparison between the self-supported and truss-supported panel concepts. In 
order for the comparison to be valid, detailed preliminary design, fabrication 
planning and test planning would be needed for several sizes in order to avoid 
overoptimistic estimates. 

Continued Investigation Is needed for film-type collectors. In 
particular, a material is needed which will satisfy the long-life, high- 
efficiency requirements of the Space Station. 
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Figure 4. Mosaic reflector of spherical segments 
assembled by RMS from Shuttle. 
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Figure 5. Assembly of integrated panel -truss modules. 




Figure 6* Sequential deployment of connected 
panel -truss modules* 
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Figure 8. Synchronously deployable concept (CREST) 
for stiff-panel reflectors. 



Figure 9. Assembly of the LPR infrared telescope.. 
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Figure 10. Deployment of preintegrated support-truss 
modules for the LPR infrared telescope. 



Figure 11 „ Support-truss for a seven-tile module of 
the LPR infrared telescope. 
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Figure 12„ Application of Pactruss to the support 
structure for the primary mirror of 
the JPL low-cost LDPx„ 


37 










85-812 



PARTIALLY DEPLOYED 


PLAN VIEW DEPLOYED 



STOWED 


Figure 13. Synchronously deployable Pactruss concept. 
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Figure 14. Deployment of conceptual model of triangular cell. 
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Elevation View 


Perspective View 


Figure 15. Views of Pactruss for offset paraboloid. 
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E ^ C eff-* 


Figure 18. Mean effective concentration ratio vs. geometric 
concentration ratio C with rms surface slope 
error a (F/D = 0.3). 



E *- C eff-* 

Figure 19. Mean effective concentration ratio vs. geometric 
concentration ratio C with rms surface slope 
error a (F/D = 0.4). 





Figure 20. Mean effective concentration ratio vs. geometric 

concentration ratio C with rms surface slope 

error a (F/D = 0.5). 
a 
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Figure 21. Mean effective concentration ratio vs. geometric 
concentration ratio C with rms surface slope 
error a (F/D = 0.6). 
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Figure 22. Mean effective concentration ratio vs. geometric 
concentration ratio C with rms surface slope 
error a (F/D = 0.7). 
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Figure 23. Effect of F/D on required geometric concentration 
ratio C with rms surface slope error a 
(E [C eff ] = 1,000). 
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Figure 26. Nominally square shallow spherical 
panel with side length "a". 
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Figure 27. Location and orientation geometry of square spherical panels. 
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Figure 28. Effect of panel orientation angle tp on maximum 
panel size a^/F vs. panel location angle e c . 
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Figure 29. Effects of pivoting the panel so that p and c 
are not coincident, on the maximum allowable 
panel size for the case where C « = 2,000, 
p = F, and \p = 45 degrees. 
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Figure 3Q. Maximum panel size a^/F vs. panel location angle 

e with different panel distances p n /F. 
c p 
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Figure 31. Maximum panel size a m /F vs. panel location angle 
e with different panel distances p_/F. 






Figure 33. PACTRUSS-supported, identical -pa 
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APPENDIX 


/* SOLARPAC.C - Program for generating the design geometry of a solar 

concentrator made up of identical square panels supported 
by a Pactruss structure. 

jmh 1/3/86 */ 


finclude "math.h" 

struct geom { 

double focal; 
double size; 
double tilt; 

double corner_coord[4] [3] ; 

} ; 

main() 

{ 

in t i , j , n , m , dta f ile,flag,N, Nnodes , Ncorners , Nmembers , Nt ota 1 ; 
int Ntypes.chan; 
char chr; 

double a,l,H,h,F, tempi, temp2,r, theta, phi, xc,yc,zc,xl,yl,zl; 

doub le *nodes , *corners , *ptemp , *pt r , *pcorn ; 

struct geom *tile_ptr; 

static double accuracy = 1.0E-12; 


a = 1.95; 

/* Panel size */ 


1 = 2.1; 

/* Bay length */ 


H = 1; 

/* Truss height */ 


h = 0; 

/* Vertical panel offset 

*/ 

F = 21.; 

/* Focal length */ 


N = 5; 

/* Radial number of panels 

*/ 


while(chr != 3) { 

printf ("\nDesign values are:\na □ %f f 1 = %f, F = %f, N = %d\n", 

a,l,F,N); 

printf(”Any changes? <none>"); 
if((chr = getcharO) == ’\r’) 
break; 

printf(”\nEnter new values separated by commas."); 
scan f ( "*F , %F , %F , !td\n" , &a , &1 , &F , &N) ; 

} 

Nnodes = (N + 1)*(N + 1); 

Ncorners = 4*N*(N + 1); 

nodes = calloc(Nnodes*3,sizeof (xc) ) ; 
corners = calloc(Ncorners*3,sizeof (xc) ) ; 
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ptemp = nodes; 

*ptemp++ = a/2 . ; 

*ptemp++ = a/2.; 

*ptemp++ = 0; 

for(i=l; i<=N; i++) { 

ptemp = nodes + 3*(i - 1); 
xl = *ptemp++; 
yl = *ptemp++; 
zl = *ptemp++; 

phi = 1 . ; 
theta = 0; 

while(fabs( theta - phi) > accuracy) { 
phi = theta; 

tempi = xl*sin(phi) - zl*cos(phi); 

r = xl*cos(phi) + zl*sin(phi) + sqrt(l*l - templ*templ) ; 
theta = asin((r - a/2.)/2./F); 

} 

*ptemp++ = r*cos ( theta) ; 

*ptemp++ = yl; 

*ptemp++ = rtsin(theta) ; 

} 

for( j=l; j<=N; j++) { 

ptemp = nodes + 3*(N + l)*j; 
yl = *(nodes + 3*j); 
zl = *(nodes + 3*j + 2); 

for(i=0; i<=N; i++) { 

*ptemp++ = *(nodes +3*i); 

*ptemp++ = yl; 

*ptemp++ = *( nodes + 3*i + 2) + zl; 

} 

} 

tile_ptr->focal = F; 
tile_ptr->size = a; 
tile_ptr->tilt = theta/2. ; 

pcorn = corners; 

for(j=l; j<=N; j++) { 

for(i=0; i<=N; i++) { 

ptr = tile_ptr->corner_coord; 
ptemp = nodes + 3*(N + l)*j + 3*i; 

*ptr++ = *ptemp++; 

*ptr++ = *ptemp++; 

*ptr — = *ptemp + h; 

Ptr — ; 

tile(tile ptr) ; 
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for(m=0; m<12; m++) 

*pcorn++ = *ptr++; 

} 

} 

Ntotal = 4*2*Nnodes +4*Ncorners + 8; 

Nmembers = 4*2*2*(N*(N + 1) - 6); 
Nmembers += 4*((N + 1)*(N + 1) - 6); 
Nmembers += 0; 

Nmembers += 4*2*(N + 1); 

Nmembers += 4*4*(N*N - 6) + 2*4*N; 
Nmembers += 12; 

Ntypes = 11; 

chan = f open ( "SOLARPAC . DTA" , "w” ) ; 


/* First the longerons */ 

/* Then the verticals */ 

/* Then the diagonals */ 

/* Then the connectors */ 

/* Then the panels */ 

/* Then the receiver; */ 

/* Separating upper and lower */ 


fprintf (chan, "^d,^d,%d\n M , Ntotal, Nmembers, Ntypes) ; 


wrtnodes ( chan , nodes , H , N) ; 
wr t_corners ( chan , corners , N) ; 
wrt_receiver(chan,F + h,a); 
wrt_members(chan,N) ; 

fclose(chan) ; 

} 


/* TILE.C - Determines the locations of the corners of a square reflector 
panel after rotation from the horizontal to the desired 
orientation. 

jmh 12/29/85 */ 

tile(ptr) 
struct geom *ptr; 

{ 

int i, j; 

double xc , yc , zc , r , tempi , temp2 , phiO , phi 1 , a , F ; 
double *pc = ptr->corner_coord; 

static double accuracy = 1.0E-10; 

a = ptr->size; 

F = ptr->focal; 
xc = *pc++ - a/2 . ; 
yc = *pc++ - a/2. ; 
zc = *pc++; 

r = sqrt(xc*xc + yc*yc) ; 

tempi = 2.*(F - zc)*r/a/(xc + yc) ; 

temp2 = 2.*r*r/a/(xc + yc) ; 


for(i=0,phil=0,phi0=l. ; (fabs (phi l-phi0)> accuracy) && (i<10000); i++) { 

53 


phiO = phil; 

phil = 0.5*atan2(temp2 + 1. - cos ( phiO ), tempi + sin(phiO)); 

} 

ptr->tilt = phil; 

tempi = a*(l. - cos(phil) )/r/r; 
temp2 = a*sin(phil)/r; 

*pc++ o xc + a/2. + xc*yc*templ; 

*pc++ = yc - a/2. + yc*yc*templ; 

*pc++ = zc - yc*temp2; 

*pc++ = xc - a/2. + xc*xc*templ; 

*pc++ = yc + a/2. + xc*yc*templ; 

*pc++ = zc - xc*temp2; 

*pc++ = xc - a/2. + xc*(xc + yc)*templ; 

*pc++ = yc - a/2. + yc*(xc + yc)*templ; 

*pc++ = zc - (xc + yc)*temp2; 

} 

wr t_nodes ( chan , nodes , H , N) 
int chan,N; 
double *nodes,H; 

{ 

int i, j, surf , quad; 

double x, y, z,*ptr,x_sign, ysign, z_delta; 

for(surf=0; surf<2; surf++) { 
zdelta = surf*H; 

for(quad=0; quad<4; quad++) ( 
ptr = nodes; 
x sign = 1. ; 
ysign =1.; 
if(quad==l ! ! quad==2) 
x_sign = -1. ; 
if(quad==2 ! ! quad==3) 
y_sign = -1. ; 

for(i=0; i<=N; i++) { 

for( j=0; j<=N; j++) { 
x = x_sign*(*ptr++) ; 
y = y_sign*(*ptr++) ; 
z = *ptr++ - zdelta; 

fprintf (chan, 10e,%. 10e,%. 10e\n" ,x,y, z) ; 


wrt_corners ( chan , corners , N ) 
int chan.N; 
double ^corners; 

{ 
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int i,j,n,quad; 

double x,y,z,x_sign,y_sign,*ptr; 

for(quad=0; quad<4; quad++) { 
ptr = corners; 
x_sign = 1 . ; 
ysign = 1. ; 
if(quad==l ! ! quad==2) 
xsign = -1. ; 
if (quad==2 ! ! quad==3) 
y_sign = -1. ; 

for( j = 1; j<=N; j++) { 
for(i=0; i<=N; i++) { 
for(n=0; n<4; n++) { 
x = x_sign*(*ptr++) > 
y = y_sign*(*ptr++) ; 
z = *ptr++; 


} 


} 


} 


fprintf(chan,"%.10e,%.10e,%. 10e\n",x,y,z) ; 

} 


wrt_receiver(chan,F,a) /* Locate corners of an axaxa cube at F */ 

int chan; 
double F,a; 

{ 

int i, j,k; 


} 


for(k=0; k<2; k++) 

for(i=0; i <2 ; i++) 

for( j=0; j <2 ; j++) 

fprintf (chan, 10e,%. 10e,%. 10e\n", (l-2*i)*a/2. , (l-2*j)*a/2. , 

F + k*a) ; 


wrt_members ( chan , N ) 
int chan,N; 

{ 

int i , j , k , n , quad , half , surf , endl , end2 , Nnodes , Ncorners , type , Cbase ; 

Nnodes = (N + 1)*(N + 1); 

Ncorners = 4*N*(N + 1); 

Cbase = 4*2*Nnodes; 

for(surf=0; surf<2; surf++) { 

for(quad=0; quad<4; quad++) { 
for(i=0; i<=N; i++) { 

for( j=l; ( j<=2*N-3-i) && j<=N; j++) { 

endl =i+(N+l)*j+ (quad + 4*surf )*Nnodes; 
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end2 = endl - N -1; /* y-wise longerons */ 

fprintf (chan, "%d,%d,%d\n", endl,end2, surf) ; 


endl = j + (N + l)*i + (quad + 4*surf )*Nnodes; 

end2 = endl - 1; /* x-wise longerons */ 

fprintf (chan, "%d,%d,%d\n" , endl, end2, surf) ; 

} 

} 

} 

for(half=0; half <2; half++) { 
for(i=0; i<=N; i++) { 

endl = i*(N + 1) + (4*surf + 2*half)*Nnodes; 

end2 = endl + Nnodes; /* x-wise connectors */ 

fprintf(chan,"%d,%d,%d\n",endl,end2,6 + surf); 

endl = i + (4*surf + 2*half)*Nnodes; 

end2 = endl - Nnodes; /* y-wise connectors */ 

if(half == 0) 

end2 = endl + 3*Nnodes; 

fprintf(chan,"%d,&d,%d\n",endl,end2,6 + surf); 

> 

} 

} 

for(quad=0; quad<4; quad++) { 
for(i=0; i<=N; i++) { 

for(j=0; j<=N && ( j<=2*N-3-i) ; j++) { 
endl = i + (N + l)*j + quad*Nnodes; 

end2 = endl + 4*Nnodes; /* Verticals */ 

fprintf (chan, "%d,%d,%d\n",endl,end2,2) ; 

} 

} 

for(j=l; j<=N; j++) { /* Panels */ 

for(i=0; i<=N && (i<=2*N-j-3) ; i++) { 
if(i>0 !! (quad==0 !! quad==2)) { 

endl = Chase + 4*i + 4*(j - 1)*(N + 1) + quad*Ncorners ; 
end2 = endl +1; 

fprintf (chan, "%d,*d, %d\n" , endl , end2 , 8) ; 
end2 = endl + 2; 

fprintf (chan, "%d,3>d,%d\n",endl,end2,8) ; 

endl += 3; 
end2 = endl - 1; 

fprintf (chan, "%d,%d,%d\n",endl,end2,8) ; 


end2 = endl - 2; 

fprintf ( chan , "%d , %d , %d\n" , endl , end2 , 8 ) ; 

} 

} 

} 

} 

for(i-0; i<4; i++) { /* Receiver */ 

endl = 8*Nnodes + 4*Ncorners + i; 
end2 = endl +4; 
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fprintf (chan, "S»d,%d,%d\n" , endl, end2, 9) ; 

} 

for(i=0; i<2; i++) { 

endl =8*Nnodes + 4*Ncorners +4*i; 
end2 = endl + 1; 

fprintf (chan, "%d,%d,%d\n", endl, end2, 9+i) ; 
end2 = endl + 2; 

fprintf (chan, "%d,%d,%d\n",endl, end2,9+i) ; 

endl += 3; 

end2 = endl - 1; 

fprintf (chan, "%d,3>d,&d\n",endl, end2, 9+i) ; 
end2 = endl - 2; 

fprintf (chan, , '%d,%d,%d\n ,, , endl, end2, 9+i) ; 

> 

} 
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